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O Abstract. We derive a Lagrangian based approach to study the compatible Hamiltonian 
structure of the dispersionless KdV and supcrsymmetric KdV hierarchies and claim that 
CNl ■ our treatment of the problem serves as a very useful supplement of the so-called r-matrix 

method. We suggest specific ways to construct results for conserved densities and Hamil- 
tonian operators. The Lagrangian formulation, via Noether's theorem, provides a method 
' to make the relation between symmetries and conserved quantities more precise. We have 

c~j . exploited this fact to study the variational symmetries of the dispersionless KdV equation. 
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t£> ■ 1 Introduction 

O 

^sO ' The equation of Korteweg and de Vries or the so-called KdV equation 

i> 
o 



in the dispersionless limit pQ 



d d d d 

■ "57 ^ e ^7 and 7^ *• e 7T" Wlth e ~^ 

at at ox ox 

reduces to 

Ut = \uu x . (1.1) 

Equation (|l.lj) . often called the Riemann equation, serves as a prototypical nonlinear partial 
differential equation for the realization of many phenomena exhibited by hyperbolic systems [2]. 
This might be one of the reasons why, during the last decade, a number of works [3] was envisaged 
to study the properties of dispersionless KdV and other related equations with special emphasis 
on their Lax representation and Hamiltonian structure. 

The complete integrability of the KdV equation yields the existence of an infinite family of 
conserved functions or Hamiltonian densities 7Y n 's that are in involution. All 7Y n 's that generate 
flows which commute with the KdV flow give rise to the KdV hierarchy. The equations of the 
hierarchy can be constructed using [3] 

u t = A n u x (x,t), n = 0,1,2,... (1.2) 
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with the recursion operator 

k = \d 2 x + u + \u x d- 1 . 

In the dispersionless limit the recursion operator becomes 

K = u+\u x d~ 1 . (1.3) 



According to (jl.2|) . the pseudo-differential operator A in (|1.3p defines a dispersionless KdV 
hierarchy. The first few members of the hierarchy are given by 



n = 

n = 1 

n = 2 

n = 3 

n = 4 



= u x , (1.4a) 

tit = §im x , (l-4b) 

u t = fu 2 u x , (1.4c) 

u t = fgtiV, (l-4d) 

u * = fi" 4 ^- ( L4e ) 

Thus the equations in the dispersionless hierarchy can be written in the general form 

u t = A n u n u x , (1.5) 

where the values of A n should be computed using (|1.3|) in (jl.2p . We can also generate A%, A2, 
A3 etc recursively using 

A n = { 1 + ^l) A n~l, n= 1,2,3,... and A = 1. 

The Hamiltonian structure of the dispersionless KdV hierarchy is often studied by taking 
recourse to the use of Lax operators expressed in the semi-classical limit [5]. In this work we 
shall follow a different viewpoint to derive Hamiltonian structure of the equations in (11. 5h . We 
shall construct an expression for the Lagrangian density and use the time-honoured method of 
classical mechanics to rederive and reexamine the corresponding canonical formulation. A single 
evolution equation is never the Euler-Lagrange equation of a variational problem. One common 
trick to put a single evolution equation into a variational form is to replace u by a potential 
function u = —w x . In terms of w, (jl.5p will become an Euler-Lagrange equation. We can, 
however, couple a nonlinear evolution equation with an associated one and derive the action 
principle. This allows one to write the Lagrangian density in terms of the original field variables 
rather than the w's, often called the Casimir potential. In Section 2 we adapt both these 
approaches to obtain the Lagrangian and Hamiltonian densities of the Riemann type equations. 
In Section 3 we study the bi-Hamiltonian structure [6]. One of the added advantage of the 
Lagrangian description is that it allows one to establish, via Noether's theorem, the relationship 
between variational symmetries and associated conservation laws. The concept of variational 
symmetry results from the application of group methods in the calculus of variations. Here one 
deals with the symmetry group of an action functional A[u] = Jq o C (x,u( n )) dx with C, the 
so-called Lagrangian density of the field u(x). The groups considered will be local groups of 
transformations acting on an open subset A4 C X U C X xU. The symbols X and U denote 
the space of independent and dependent variables respectively. We devote Section 4 to study 
this classical problem. Finally, in Section 5 we make some concluding remarks. 

2 Lagrangian and Hamiltonian densities 

For u = —w x (|1.5p becomes 

w xt = A n (-l) n w2w 2x . (2.1) 
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The Frechet derivative of the right side of (|2.ip is self-adjoint. Thus we can use the homotopy 
formula [7] to obtain the Lagrangian density in the form 



A i'_iy i + 1 



In writing (|2.2p we have subtracted a gauge term which is harmless at the classical level. The 
subscript n of C merely indicates that it is the Lagrangian density for the nth member of the 
dispersionless KdV hierarchy. The corresponding canonical Hamiltonian densities obtained by 
the use of Legendre map are given by 

Hn = ^ -u n+2 . (2.3) 

(n + l)(n + 2) V ; 

Equation (|1.5[) can be written in the form 



u, + ^l = ,,4) 

with 

There exists a prolongation of (j 1 . 5 j) or (|2.4p into another equation 

v t+ 6 -^l = 0, v = v(x,t) (2.6) 
ou 

with the variational derivative 

& v^, d k d d k u 



OU OX K OUh x ox^ 

such that the coupled system of equations follows from the action principle [8] 
/ C c dxdt = 0. 



The Lagrangian density for the coupled equations in (|2.4p and (|2.6p is given by 

C c = \{vu t - uv t ) - p[u]v x . 
For p[u] in ([53]) . (f2"U|) reads 

v t = A n u n v x . (2.7) 
For the system represented by (|1.5p and (|2.7p we have 



C c n = \(vut - uvt) + -^-u n+1 v x . (2.8) 

The result in (|2.7p could also be obtained using the method of Kaup and Malomed [S]. Referring 
back to the supersymmetric KdV equation [10] we identify f as a fermionic variable associated 
with the bosonic equation in (jl.5p . It is of interest to note that the supersymmetric system is 
complete in the sense of variational principle while neither of the partners is. The Hamiltonian 
density obtained from the Lagrangian in (12. 8p is given by 

It remains an interesting curiosity to demonstrate that the results in (|2.3p and (|2.9p represent 
the conserved densities of the dispersionless KdV and supersymmetric KdV flows. We demon- 
strate this by examinning the appropriate bi-Hamiltonian structures of (|1.5|) and the pair (|1.5|) 
and (Ej 
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3 Bi-Hamiltonian structure 

Zakharov and Faddeev [11] developed the Hamiltonian approach to integrability of nonlinear 
evolution equations in one spatial and one temporal (1+1) dimensions and Gardner [12], in 
particular, interpreted the KdV equation as a completely integrable Hamiltonian system with d x 
as the relevant Hamiltonian operator. A significant development in the Hamiltonian theory is 
due to Magri [6] who realized that integrable Hamiltonian systems have an additional structure. 
They are bi-Hamiltonian, i.e., they are Hamiltonian with respect to two different compatible 
Hamiltonian operators. A similar consideration will also hold good for the dispersionless KdV 
equations and we have 



= d x ( = \ (ud x + d x u) ( 1 ) , n = 1, 2, 3 . . . . 



Here 

H = I Udx. (3.2) 



It is easy to verify that for n = 1, (|2.3p . (|3.ip and (|3.2p give (jl.4b|) . The other equations of the 
hierarchy can be obtained for n = 2,3,4, ... . The operators T>\ = d x and T> 2 = \ (ud x + d x u) 
in (|3.ip are skew-adjoint and satisfy the Jacobi identity. The dispersionless KdV equation, in 
particular, can be written in the Hamiltonian form as 

u t = {u(x),H 1 } 1 and u t = {u(x), H } 2 

endowed with the Poisson structures 

{u{x),u(y)}x=V 1 8(x -y) and {u{x), u(y)} 2 = V 2 5(x - y). 

Thus T>\ and T> 2 constitute two compatible Hamiltonian operators such that the equations 
obtained from (II. 5p are integrable in Liouville's sense [6]. Thus H n 's in (|2.3p via (|3.2p give 
the conserved densities of (jl.5p . In other words, W n 's generate flows which commute with the 
dispersionless KdV flow and give rise to an appropriate hierarchy. It will be quite interesting 
to examine if a similar analysis could also be carried out for the supersymmetric dispersionless 
KdV equations. 

The pair of supersymmetric equations ut = u n u x and vt = u n v x can be written as 



where rj = ( U ^ , = ^ and H% = f H c n dx. In (|3.3p Ji and J2 stand for the matrices 

and ,,= C °. (3.4) 



-1 J z \ -u / 

Since for different values of n represent the conserved Hamiltonian densities obtained by the 
use of action principle, the supersymmetric dispersionless KdV equations will be bi-Hamiltonian 
provided Ji and J 2 constitute a pair of compatible Hamiltonian operators. Clearly, Ji and J 2 
are skew-adjoint. Thus Ji and J 2 will be Hamiltonian operators provided we can show that [5] 



Prvj^ejj = 0, * = 1,2. 



(3.5) 
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Here pr stands for the prolongation of the evolutionary vector field v of the characteristic Ji9. 
The quantity prvj^ is calculated by using 

Prv j4 , = X;^(E( J «Vn^. D i = lL> M,"=l,2. (3.6) 



drf- cte? 



In our case the column matrix 9 = I ^ J represents the basis univectors associated with the 

variables 77 = ^ ^ J . Understandably, 9 U and ?] M denote the components of 9 and r/ and (Jj)^ 
carries a similar meaning. The functional bivectors corresponding to the operators Jj is given by 

Sj, = \J 9 T AJi9dx (3.7) 

with 9 T , the transpose of 9. From (|3.4p . ()3.6p and (|3.7p we found that both Ji and J 2 satisfy ()3.5p 
such that each of them constitutes a Hamiltonian operator. Further, one can check that Ji and J 2 
satisfy the compatibility condition 

P r v Jl0 (e j2 ) + prv j20 (e Jl ) = 0. 



This shows that (13. 3f) gives the bi-Hamiltonian form of supersymmetric dispersionless KdV 
equations. The recursion operator defined by 



reproduces the hierarchy of supersymmetric dispersionless KdV equation according to 
rj t = A n r) x . 

for n = 0, 1, 2, . . . . This verifies that -r^'s as conserved densities generate flows which commute 

sin 

with the supersymmetric dispersionless KdV flow. 



4 Variational symmetries 

The Lagrangian and Hamiltonian formulations of dynamical systems give a way to make the re- 
lation between symmetries and conserved quantities more precise and thereby provide a method 
to derive expressions for the conserved quantities from the symmetry transformations. In its 
general form this is referred to as Noether's theorem. More precisely, this theorem asserts that 
if a given system of differential equations follows from the variational principle, then a continu- 
ous symmetry transformation (point, contact or higher order) that leaves the action functional 
invariant to within a divergence yields a conservation law. The proof of this theorem requires 
some knowledge of differential forms, Lie derivatives and pull-back [5]. We shall, however, carry 
out the symmetry analysis for the dispersionless KdV equation using a relatively simpler mathe- 
matical framework as compared to that of the algebro- geometric theories. In fact, we shall make 
use of some point transformations that depend on time and spatial coordinates. The approach 
to be followed by us has an old root in the classical-mechanics literature. For example, as early 
as 1951, Hill [13] provided a simplified account of Noether's theorem by considering infinitesimal 
transformations of the dependent and independent variables characterizing the classical field. 
We shall first present our general scheme for symmetry analysis and then study the variational 
or Noether's symmetries of the dispersionless KdV equation. 
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Consider the infinitesimal transformations 

x*' = x i + 5x\ 5x i = eC(x, f) (4.1a) 

and 

f' = f + Sf, 6f = er ) (xJ) (4.1b) 

for a field variable / = f(x,t) with e, an arbitrary small quantity. Here x = {x ,^ 1 }, x° = t 
and x 1 = x. Understandably, our treatment for the symmetry analysis will be applicable to 
(1 + 1) dimensional cases. However, the result to be presented here can easily be generalized to 
deal with (3 + 1) dimensional problems. For an arbitrary analytic function g = g(x l ,f), it is 
straightforward to show that 

dg = tXg 

with 

the generator of the infinitesimal transformations in (|4.ip . A similar consideration when applied 
to h = h{x\ f, fi) with fi = J£ gives 

5h = eX'h (4.3) 

with 

jr = X (4.4) 

Understandably, X' stands for the first prolongation of X. To arrive at the statement for the 
Noether's theorem we consider among the general set of transformations in (|4.ip only those that 
leave the field-theoretic action invariant. We thus write 

L{x\ /, fi)d(x) = £{x l \ /', //MA (4-5) 

where d(x) = dxdt. In order to satisfy the condition in (|4,5p we allow the Lagrangian density 
to change its functional form C to £ . If the equations of motion, expressed in terms of the 
new variables, are to be of precisely the same functional form as in the old variables, the two 
density functions must be related by a divergence transformation. We thus express the relation 
between £' and C by introducing a gauge function B % (x, f) such that 



C'ix 1 ', /', fi')d^) = Ctf', /', fi')d(x') - e^—jd(x') + o{e 2 ). (4.6) 



dx l 

The general form of (|4.6p for the definition of symmetry transformations will allow the scale 
and divergence transformations to be considered as symmetry transformations. Understandably, 
the scale transformations give rise to Noether's symmetries while the scale transformations in 
conjunction with the divergence term lead to Noether's divergence symmetries. Traditionally, 
the concept of divergence symmetries and concommitant conservation laws are introduced by 
replacing Noether's infinitesimal criterion for invariance by a divergence condition p3]. However, 
one can directly work with the conserved densities that follow from (|4.6p because nature of the 
vector fields will determine the contributions of the gauge term. For some of the vector fields 
the contributions of B l to conserved quantities will be equal to zero. These vector fields are 
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Noether's symmetries else we have Noether's divergence symmetries. In view of (|4.5p . (|4.6p can 
be written in the form 



C(x l \ f, h')d{x') = C(x\ f, fr)d(x) + e^-dix). 
Again using C for h in (|4.3p . we have 

C{j', /', = /, /O [d(x) + edf (re, ft)] + eX'C{x\ f, fi)d(x) 

From (|4,7p and (|4.8p . we write 

dx % dx l 

Using the value of X' from (|4.4p in (|4.9p , is obtained in the final form 
dxi - dx* + * dx* + V df + ^ * ^ 9/, • 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



Thus we find that the action is invariant under those transformations whose constituents £ and r/ 
satisfy ()4.10|) . The terms in (I4.10|) can be rearranged to write 



d 



dC 



— \B*-ec+ (?f s -r,) — \ + {eh - v) 



dfi 





d 


(-)] 




dx l 





0. 



(4.11) 



The expression inside the squared bracket stands for the Euler-Lagrange equation for the clas- 
sical field under consideration. In view of this, (|4.1ip leads to the conservation law 



dV 

dx % 







with the conserved density given by 

r = B>-e£+(Pf j - v )—. 



(4.12) 



(4.13) 



In the case of two independent variables (x ,^ 1 ) = (t, x), (|4.12p can be written in the explicit 
form 



dl° dl 1 
+ =0. 

dt dx 



(4.14) 



From (|2.2p the Lagrangian density for the dispersionless KdV equation is obtained as 

C = \w t w x + \wl. (4.15) 
Identifying / with w we can combine (|4.13p . (|4.14p and (14.151) to get 



B? + w t Bl - \$wl - j^w t wl + \i}w 2 x + \iiw t w 2 x - \r} t w x - r] w w t w x 

+ Bl + w x Bl + \£wl + \iiw x + \i° x w 2 t + \w x wUi + ltlw 2 x w t - \j] x wl 
- hvxWt - hwwl + jCwlw t = 0. 



(4.16) 
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In writing (|4.16p we have made use of (|2.ip with n = 1. Equation (|4.16p can be globally satisfied 
iff the coefficients of the following terms vanish separately 



w x or 


w° t : B t ° + B\ = 0, 


(4.17a) 


Wt '■ 


R° - \j] x = 0, 


(4.17b) 


w 2 ■ 

t 




(4 17c) 


^x 


- k = o, 


(4.17d) 


w 2 , : 




(4.17e) 


ufi • 


4St 4 'm> ^ 2^£ ' 


(4 17f) 


w i ■■ 


if 1 - 


(4.17g) 


w t w x : 


= o, 


(4.17h) 


w t w 2 x : 


2^ui t 4S2 u i 


(4.17i) 


w t wl : 


±£° = 


(4-17j) 


VJt w x ■ 


±£° = 


(4.17k) 



Equations in (|4.17p will lead to finite number of symmetries. This number appears to be 
disappointingly small since we have a dispersionless KdV hierarchy given in (jl.5p . Further, 
symmetry properties reflecting the existence of infinitely many conservation laws will require 
an appropriate development for the theory of generalized symmetries. In this work, however, 
we shall be concerned with variational symmetries only. 

From ffiTcD . gTfl and KTfB we see that £° is only a function of t. We, therefore, write 

e(x,t,w)=(3(t). (4.18) 

Also from ( 4,17gP , (|4.17ip and (|4,18p we see that £ is not a function of w. In view of (|4.17hp 
and (HHHD, flUSH) gives 

£* - |A = 
which can be solved to get 

Z 1 = lP t x + a(t), (4.19) 
where a(t) is a constant of integration. Using (|4.19|) in (|4.17ep we have 

rjx = \Pttx + \a t . (4.20) 
The solution of (|4.20p is given by 

V = iPttx 2 + \a t x + 7 (t) (4.21) 
with 7(t), a constant of integration. In view of (|4.2ip . (|4.17bp and (|4.17dp yield 

B° = lp tt xw + \a t w (4.22) 

and 

B 1 = ^f3 tt tx 2 w + \a tt xw. (4.23) 
Equations (I4.22|) and (|4.23p can be combined with H4. 17aj) to get finally 

Pttt = and a tt = 0. (4.24) 
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From (|4,24p we write 

= \ait 2 + a 2 t + a 3 (4.25) 

and 

a = b 1 t + b 2 , (4.26) 

where a's and Vs are arbitrary constants. Substituting the values of (5 and a in (|4.18p . (|4.19|) . 
(|4.2ip we obtain the infinitesimal transformation, £°, £} and rj, as 

f = l^t 2 + a 2 t + a 3 , (4.27a) 

t 1 = \ {ait + a 2 )x + bit + b 2 , (4.27b) 

= \aix 2 + §bix + b 3 . (4.27c) 

In writing (|4.27cp we have treated j(t) as a constant and replaced it by 63. Implication of this 
choice will be made clear while considering the symmetry algebra. In terms of (|4.27p . (|4.2p 
becomes 

X = aiVi + a 2 V 2 + a 3 V 3 + 61V4 + b 2 V 5 + 63^, 

where 



F 3 = |, F 4 = i|- + §xA Fs = |_, % * (4.28) 
at ax J au> ax aiu 

It is easy to check that the vector fields V\,...,Vq satisfy the closure property. The commutation 
relations between these vector fields are given in Table [TJ 



Table 1. Commutation relations for the generators in (|4.28[) . Each element Vij in the Table is represented 
byV ij = \y i ,V j \. 







Vi 


v 3 


v A 


v 5 


v 6 







-Vi 


-v 2 










v 2 


Vi 





-v 3 




-\Vs 





v 3 


v 2 


^3 





v 5 








v A 





-\v A 


-v 5 





-1^6 





v 5 


\v, 















v 6 





















The symmetries in (I4.28P are expressed in terms of the velocity field and depend explicitly 
on x and t. Looking from this point of view the symmetry vectors obtained by us bear some 
similarity with the so called 'addition symmetries' suggested independently by Chen, Lee and 
Lin [15] and by Orlov and Shulman [16] . It is easy to see that V 2 to V% correspond to scaling, 
time translation, Galilean boost, space translation and translation in velocity space respectively. 
The vector field V\ does not admit such a simple physical realization. However, we can write V\ 
as Vi = \tV 2 + \xV 4 . 

Making use of P~T5]) . (|£2Zj) . flj]25D , (|4T23|) and (|4T2"rjj) we can write the expressions for the 
conserved quantities in (|4.13j) as 

1° = \aixw + \b x w - \i°wl + \^w 2 x - \r]w x , (4.29a) 
X 1 = ±£V + l^wtwl + ±ewl - \rm ~ h w l ( 4 - 29b ) 
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The expressions for 2° and I 1 are characterized by £ l and 77, the values of which change as we 
go from one vector field to the other. The first two terms in X° stand for the contribution of B° 
and there is no contribution of the gauge term in X 1 since from (|P31) and (Oij) B 1 = 0. For 
a particular vector field a\ and b\ may either be zero or non zero. One can verify that except for 
vector fields V\ and V4, a\ = b\ = such that V2, V3, V5 and Vq are simple Noether's symmetries 
while V\ and V4 are Noether's divergence symmetries. Coming down to details we have found 
the following conserved quantities from (|4.29ap and (|4.29bl) 





= — it 2 W? + \xtvJt — -r^X^Wrr, 
6 8 x k x 12 x ' 


(4.30a) 


T 1 

VI 


1 x 3 1 3j.2 2 1 lj.2 2 12 122 
= -gXtW x + WtW x + ^1 — yrjX — Wj,, 


(4.30b) 


T° 

v 2 


— —I^ 3 _i_ \xw 2 . 


(4.30c) 


T 1 




(4.30d) 




= -±u; 3 

4 '"1! 


(4.30e) 


T 1 




(4.30f) 




= t$w + \tw 2 x — \xw x , 


(4.30g) 


T 1 


= \tw x — \xwt — \xw 2 xl 


(4.30h) 




= ±w 2 

2 w xi 


(4.30i) 




= W 
2 xi 


(4.30j) 




— 2^xi 


(4.30k) 


T 1 


1 S 2 
= ~2^ - I W x- 


(4.301) 



It is easy to check that the results in (|4,30p is consistent with (|4,14p . The pair of conserved 
quantities corresponding to time translation, space translation and velocity space translation, 
namely, { (l4TTO . (l430f]) }. { (|430T1) ,( [430jT )i and { (|430kl) . (l430ll) l do not involve x and t explicitly. 



Each of the pair in conjunction with (|4.14p give the dispersionless KdV equation in a rather 
straightforward manner. As expected (|4.30ep stands for the Hamiltonian density or energy 
of Ob) . 



5 Conclusion 

Compatible Hamiltonian structures of the dispersionless KdV hierarchy are traditionally ob- 
tained with special attention to their Lax representation in the semiclassical limit. The deriva- 
tion involves judicious use of the so-called r-matrix method [T7]. We have shown that the 
combined Lax representation-r-matrix method can be supplemented by a Lagrangian approach 
to the problem. We found that the Hamiltonian densities corresponding to our Lagrangian rep- 
resentations stand for the conserved densities for the dispersionless KdV flow. We could easily 
construct the Hamiltonian operators from the recursion operator which generates the hierarchy. 
We have derived the bi-Hamiltonian structures for both dispersionless KdV and super symmetric 
KdV hierarchies. As an added realism of the Lagrangian approach we studied the variational 
symmetries of equation (|1.4bp . We believe that it will be quite interesting to carry out similar 
analysis for the super symmetric KdV pair in (fl~4b|) and for n = 1 limit of (f277|> . 



Acknowledgements 



This work is supported by the University Grants Commission, Government of India, through 
grant No. F.32-39/2006(SR). 



Lagrangian Approach to Dispersionless KdV Hierarchy 



11 



References 

[1] Zakharov V.E., Benney equations and quasiclassical approximation in the method of the inverse problem, 
Fund. Anal. Appl. 14 (1980), 89-98. 

[2] Olver P.J., Nutku Y., Hamiltonian structures for systems of hyperbolic conservation laws, J. Math. Phys. 
29 (1988), 1610-1619. 

Brunelli J.C., Dispersionless limit of integrable models, Braz. J. Phys. 30 (2000), 455-468, jnTm.SI/0207042| 



[3] Arik M., Neyzi F., Nutku Y., Olver P.J., Verosky J.M., Multi-Hamiltonian structure of the Born-Infeld 
equation, J. Math. Phys. 30 (1989), 1338-1344. 

Das A., Huang W.J., The Hamiltonian structures associated with a generalized Lax operator, J. Math. Phys. 
33 (1992), 2487-2497. 

Brunelli J.C., Das A., Properties of nonlocal charges in the supersymmetric two boson hierarchy, Phys. 



[5] 
[6] 

[r 

[9 
[10 

[11 

[12 

[13 

[14 
[15 

[16 

[ir 



Lett. B 354 (1995), 307-314, hep-th/9504030 



Brunelli J.C., Das A., Supersymmetric two-boson equation, its reductions and the nonstandard supersym- 
metric KP hierarchy, Intern. J. Modern Phys. A 10 (1995), 4563-4599, |hep-th/9505093| 
Brunelli J.C., Hamiltonian structures for the generalized dispersionless KdV hierarchy, Rev. Math. Phys. 8 
(1996), 1041-1054, [solv-int/960100T] 

Brunelli J.C., Das A., A Lax description for polytropic gas dynamics, Phys. Lett. A 235 (1997), 597-602, 

|solv-int/9706005| 

Brunelli J.C., Das A., The sTB-B hierarchy, Phys. Lett. B 409 (1997), 229-238, |hep-th/9704126] 



Brunelli J.C., Das A., A Lax representation for Born-Infeld equation, Phys. Lett. B 426 (1998), 57-63, 
|hep-th/9712081| 

[4] Calogero F., Degasperis A., Spectral transform and soliton, North-Holland Publising Company, New York, 
1982. 



Olver P.J., Application of Lie groups to differential equation, Springer- Verlag, New York, 1993. 

Magri F., A simple model of the integrable Hamiltonian equation, J. Math. Phys. 19 (1978), 1156-1162. 

Frankel T., The geometry of physics, Cambridge University Press, UK, 1997. 

Ali Sk.C, Talukdar B., Das U., Inverse problem of variational calculus for nonlinear evolution equations, 



Acta Phys. Polon. B 38 (2007), 1993-2002, nlin. SI/0603037 



Kaup D.J., Malomed B.A., The variational principle for nonlinear waves in dissipative systems, Phys. D 87 
(1995), 155-159. 

Barcelos-Neto J., Constandache A., Das A., Dispersionless fermionic KdV, Phys. Lett. A 268 (2000), 342- 



351, solv-int/9910001 



Zakharov V.E., Faddeev L.D., Korteweg-de Vries equation: a completely integrable Hamiltonian systems, 
Fund. Anal. Appl. 5 (1971), 18-27. 

Gardner C.S., Korteweg-de Vries equation and generalizations. IV. The Korteweg-de Vries equation as 
a Hamiltonian system, J. Math. Phys. 12 (1971), 1548-1551. 

Hill E.L., Hamilton's principle and the conservation theorems of mathematical physics, Rev. Modern Phys. 
23 (1951), 253-260. 

Gelfand I.M., Fomin S.V., Calculus of variations, Dover Publ., 2000. 

Chen H.H., Lee Y.C., Lin J.E., On a new hierarchy of symmetries for the Kadomtsev-Petviashvilli equation, 
Phys. D9 (1983), 439-445. 

Orlov A.Yu., Shulman E.I., Additional symmetries for integral and conformal algebra representation, Lett. 
Math. Phys. 12 (1986), 171-179. 

Faddeev L.D., Takhtajan L.A., Hamiltonian methods in the theory of solitons, Springer, Berlin, 1987. 



